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O ■ Abstract 

! Pion propagation in a hadronic fluid with a non-homogeneous relativistic flow is 

I studied in terms of the hnear sigma model. The wave equation turns out to be equiv- 

alent to the equation of motion for a massless scalar field propagating in a curved 
Qi spacetime geometry. The metric tensor depends locally on the soft pion dispersion re- 

^ [ lation and the four- velocity of the fluid. For a relativistic flow in curved spacetime the 

apparent and trapping horizons may be deflned in the same way as in general relativity. 
An expression for the analog surface gravity is derived. 

>: 

^ ! 1 Introduction 

00; 

^ ■ In current understanding the matter created in heavy-ion collisions behaves as a nearly 

^ . perfect expanding fluid [1] under extreme conditions of very high density and tempera- 

CN ■ ture. This hydrodynamic behavior was observed at Brookhaven's Relativistic Heavy Ion 

^ I Collider (RHIC) and recently confirmed by the ALICE collaboration in Pb-Pb collisions at 

the LHC [2, 3]. In high energy collisions the produced particles are predominantly pions. 
^ ■ The agreement between the pion production results reported in [2] and the theoretical hy- 

drodynamical model predictions [4] is truly remarkable. A realistic hydrodynamic model 
may be constructed [5] in which a transverse expansion is superimposed on a longitudinal 
boost invariant expansion [6]. 

It is often stated by particle physicists that heavy ion collisions create mini Big Bangs^ 
- events in which matter is created under extreme conditions of high density and high 
temperature resembling the conditions in the early universe a fraction of a second after the 
Big Bang. The expansion of hadronic matter that takes place immediately after a heavy 
ion collision has certain similarity with the cosmological expansion. However, the analogy 
is rather superficial since in a cosmological expansion of spacetime after the Big Bang the 



few recent quotations from the press: "What we're doing is reproducing the conditions that existed at 
the very early universe, a few millionths of a second after the Big Bang," [7] ; "The Large Hadron Collider has 
successfully created a "mini-Big Bang" by smashing together lead ions instead of protons." [8]; "The collisions 
generated mini Big Bangs and the highest temperatures and densities ever achieved in an experiment" [9]. 
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gravity plays the essential role, whereas high energy collisions and subsequent expansion 
does not involve gravity at all. Although the Minkowski spacetime with expanding hadronic 
matter can be mapped into an expanding spacetime, the resulting spacetime is still fiat. 
However, we will demonstrate here that in high energy collisions a much closer analogy 
with cosmology may be drawn owing to the effective analog gravity with essentially curved 
geometry. Various aspects of analog gravity (for a review and extensive list of references see 
[10]) have been studied in acoustics [11] optics [12], superfluidity [13], black hole accretion 
[14, 15], and hadronic fluid near the QCD chiral phase transition [16]. In this paper we 
study in detail the framework of analog gravity provided by a hadronic fluid at nonzero 
temperature for the whole range of temperatures below the chiral phase transition. We show 
that the analog cosmological spacetime corresponds to a contracting FRW universe with a 
nontrivial apparent horizon. 

Strongly interacting matter is described at the fundamental level by a nonabelian gauge 
theory called quantum chromo dynamics (QCD). At large distances or small momenta, the 
QCD exhibits the phenomena of quark confinement and chiral symmetry breaking. At low 
energies, the QCD vacuum is characterized by a nonvanishing expectation value [17]: (■0V') ~ 
(235 MeV)"', the so called quark condensate, which describes the density of quark- ant iquark 
pairs found in the QCD vacmim and its nonvanishing value is the manifestation of chiral 
symmetry breaking. The phenomenological importance of the chiral transition and possible 
experimental signatures have been discussed by Harris and Miiller [18]. 

The chiral symmetry breaking and restoration at finite temperature may be conveniently 
stTidicd using the linear sigma model [19, 20] originally proposed as a model for strong nuclear 
interactions [21]. Today, the linear sigma model serves as an effective model for the low- 
energy (low-temperature) phase of QCD. The basic model involves four scalar fields (three 
pions and a sigma meson) and two-fiavor constituent quarks. In the chirally symmetric phase 
at temperatures above the chiral transition point the mesons are massive with equal masses 
and quarks are massless. In the chirally broken phase the pions are massless, whereas the 
quarks and sigma meson acquire a nonzero mass proportional to the chiral condensate. At 
temperatures below the chiral phase transition point the pions, although being massless, 
propagate slower than light [22, 23, 24] with a velocity approaching zero at the critical 
temperature. Hence, it is very likely that there exists a region where the fiow velocity 
exceeds the pion velocity and the analog trapped region may form. In our previous paper 
[16] we have demonstrated that a region containing analog trapped surfaces forms near the 
chiral phase transition. The purpose of this paper is to study general conditions for the 
formation of a trapped region with the inner boundary as a marginally trapped surface 
which we refer to as the analog apparent horizon. Our approach is based on the linear sigma 
model combined with a boost invariant Bjorken type spherical expansion. A similar model 
has been previously studied in the context of disoriented chiral condensate [25]. 

The remainder of the paper is organized as follows. In Sec. 2 we describe the properties 
and the dynamics of the chiral fluid at flnite temperature. The analog geometry of the 
expanding chiral fluid is studied in Sec. 3 in which we derive the condition for the analog 
apparent horizon and study the analog Hawking effect. In the conchiding section. Sec. 4, we 
summarize our results and discuss physical consequences. Finally, in Appendix A we outline 
basic notions related to trapped surfaces in black hole physics and cosmology. 
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2 Chiral fluid 



In this section we focus on physics of hadrons at finite temperature and study the properties 
and the dynamics of an expanding chiral fluid. We base our study on a linear sigma model 
with no fermions which we describe in Sec. 2.1. In Sec. 2.2 we calculate the effective velocity 
of pions propagating in a chiral medium. We model the fluid expansion on a boost invariant 
spherical expansion of the Bjorken type which we describe in Sec. 2.3 

2.1 Linear sigma model 

Consider a linear sigma model at finite temperature in a general curved spacetime back- 
ground. For our purpose it is sufficient to study the model with no constituent fermions. 
The thermal bath provides a medium which may have an inhomogeneous velocity field. 
The dynamics of mesons in such a medium is described by an effective chirally symmetric 
Lagrangian of the form [26] 

1 171 ^ A 

where u^J^ is the velocity of the fluid, and (^^j^ is the background metric. The mesons </? = (cr, 
tt) constitute the (|, |) representation of the chiral SU(2) xSU(2). The parameters a and h 
depend on the local temperature T and on the parameters of the model mo and A and may 
be calculated in perturbation theory. At zero temperature the medium is absent in which 
case a — \ and 6 = 0. 

If mo < the chiral symmetry will be spontaneously broken. At the classical level, the 
a and tt fields develop nonvanishing expectation values such that at zero temperature 

H^ + (^)^ = -^-/.^- (2) 

It is convenient to choose here 

(tt,) = 0, ((7) = U- (3) 

At nonzero temperature the expectation value (cr) is temperature dependent and vanishes 
at the chiral transition point. Redefining the fields 

^ + (^'(x) = ((7, tt) + ((7'(x), 7r'(x)), (4) 

where tt' and a' are quantum fluctuations around the constant values tt = and a — (cr) 
respectively, we obtain the effective Lagrangian in which the chiral symmetry is explicitly 
broken: 

= \{ag-^ + bu'^und,^'d.^' - ^a" - ^tt" - gaV' - ^i^")'. (5) 

The fields a' and tt' correspond to the physical sigma meson and pions, respectively. The 
effective masses and the trilinear coupling g are functions of a defined as 

m^ = mo + 3A(j^, 

ml = "^0 + -^^^^ (6) 
g = Act. 
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For temperatures below the chiral transition point the meson masses are given by 



m^ = 0; mi = 2\a\ (7) 

in agreement with the Goldstone theorem. The temperature dependence of the chiral con- 
densate a is obtained by minimizing the thcrmodynamical potential Vt = ~{T/V) \nZ with 
respect to a at fixed inverse temperature /3. At one loop order, the extremum condition 
reads [20] 

'^^ = - ^ ^^^'"^^ " i ^^^'^'^^ ' 

where 

uj^^\p\; uja^ip" + mlf'^ (9) 
are the energies of the tt and a particles respectively and 

^M^) - (10) 

is the Bose-Einstein distribution function. Eq. (8) has been derived from the zero-order 
thermodynamical potential with meson masses at one loop order [20]. This approximation 
corresponds to the leading order va.l/N expansion, where N is the number of scalar fields 
[27]. In our case, iV = 4. The right-hand side of (8) depends on a through the mass m„ given 
by (7). The behavior of a near the critical temperature should be analyzed with special care. 
A straightforward solution to (8) as a function of temperature exhibits a weak first-order 
phase transition [20,28]. However, Pisarski and Wilczek have shown on general grounds that 
the phase transition in SU(2)xSU(2) chiral models should be of second order [29]. Hence, it 
is generally believed that a first-order transition is an artifact of the one loop approximation. 
Two loop calculations [30] make an improvement and confirm the general analysis of [29] . It 
is possible to mimic the second order phase transition even with (8) by making the cr-meson 
mass temperature independent all the way up to the critical temperature and equal to its 
zero-temperature mean field value given by 

ml = 2\fl (11) 

instead of (7). We fix the coupling A from the values of m,^ and /„■ for which we take m„ — 1 
GeV and /^^ = 92.4 MeV as a phenomenological input. In Fig. 1 we plot the solutions to 
(8) for both temperature dependent and temperature independent 711^ exhibitting apparent 
first and second order phase transitions, respectively. In the rest of our paper we employ the 
solution that corresponds to the second order phase transition. For our choice of parameters 
we find numerically = 182.822 MeV. 

The propagation of pions is governed by the equation of motion 



where 



(a g^"" + h u^'u'')d^'K\ + 1/((j, 7r)7r = 0, (12) 
y (a, tt) = + (7(7 + A((7' + tt") (13) 
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50 100 150 7^ 



Figure 1: Chiral condensate as a function of temperature for a temperature independent (full 
line) and temperature dependent m„ (dashed line) representing 2nd order phase transition 
and first order (discontinuous) phase transition, respectively. The critical temperature of 
the second order phase transition is indicated by T^. 



is the interaction potential. In the comoving reference frame in flat spacetime, equation (12) 
reduces to the wave equation 

(d^ - clA + ^V)7r = 0, (14) 
a 

where the quantity defined by 

+ (15) 

is the pion velocity. As we shall demonstrate in the next section. The constants a and b may 
be derived from the finite-temperature perturbation expansion of the pion self energy. 



2.2 Pion velocity 

At temperatures below the chiral transition point the pions are massless. However, the 
velocity of massless particles in a medium is not necessarily equal to the velocity of light - 
in the chiral fluid pions usually propagate slower than light^. The pion velocity in a sigma 

^If the chiral fermions are present pions become superluminal in certain range of temperature and baryon 
chemical potential [26] . 
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model at finite temperature has been calculated at one loop level by Pisarski and Tytgat in 
the low temperature approximation [22] and by Son and Stephanov for temperatures close 
the chiral transition point [23, 24]. It has been found that the pion velocity vanishes as one 
approaches the critical temperature. Here we summarize the calculation of the parameters 
a and b in the entire range of temperatures in the chiral symmetry broken phase [26] . 

The pion velocity may be derived from the self energy E(g, T) in the limit when the 
external momentum q approaches 0. For a flat background geometry g^^^ = 7)^^, the inverse 
pion propagator derived directly from the effective Lagrangian (5) as 



may in the limit g ^ be expressed in the form 



Z,A-^ = q^q, - -q^q^ 



d d 



[S(g,T)-E(g,0)) 



+ 



(16) 



(17) 



J?=o 



where the ellipsis denotes the terms of higher order in g'^. The g'^ independent term of the self 
energy absorbs in the renormalized pion mass, equal to zero in the chiral symmetry broken 
phase. The subtracted T = term has been absorbed in the wave function renormalization 
factor Zt,. By comparing this equation with Eq. (16) written in the comoving frame as 



A-' = (a + h)qi - aq^ 



(18) 



we can express the parameters a and 6, and hence the pion velocity, in terms of second 
derivatives of S(g,T) evaluated at = 0. At one loop level the only diagram that gives a 
nontrivial q-dependence of S is the bubble diagram. Subtracting the T = term one finds 
[24] 
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E(g,T)-E(g,0) = -4//A^ 
|[ns(a;^) + nB{uJa,q)] 



1 3 2u;T,2u;a^q 



1 



+ — 

^ 1 1 



+ 



90, 



(19) 



where cUg-.q = [(p — 9)^ + wi^]^^^. Here we take m„ to be a function of a through Eq. (7). A 
straightforward evaluation of the second derivatives of S(g) at g^ = yields 



a = 1 + 



165^^ /■ d^p 



rrvi 



3 rni 



'nB{oJn) riBiuJa) 



ml, 



I 



(27r)3 _ 
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(20) 



(21) 



The pion velocity c,r ^ given by (15) depends on temperature explicitly through the thermal 
distribution function ub and implicitly through the chiral condensate a given by Eq. (8). In 
Fig. 2 we plot Cjr as a function of temperature corresponding to two solutions depicted in 
the left panel. 
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I -/^[MeV] 

50 100 150 Tc 

Figure 2: Pion velocity as a function of temperature for a temperature independent (full 
line) and temperature dependent m^r (dashed line). The critical temperature of the second 
order phase transition is indicated by Tc. 



2.3 Spherical Bjorken expansion 

In order to explore the analogy between the chiral-fluid and cosmological expansions we 
consider a boost invariant spherically symmetric Bjorken type expansion [6] in Minkowski 
background spacetime. In radial coordinates = [t, r, -t}, ip) the fluid velocity is given by 

n^ = (7,7^,0,0) = (t/r,r/r,0,0), (22) 

where v = r/t is the radial three- velocity and r = y/t"^ — r"^ is the proper time. Using the so 
called radial rapidity 

1 t+r , , 

the velocity is expressed as 

= (cosh y, sinh y,0,0), (24) 
and hence, the radial three-velocity is 

V = tanhy. (25) 
It is convenient to change {t, r, i), ip) to new coordinates (r, y, ^, ip) via the transformation 

if: = r coshy, 

r = T sinh y. (26) 
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In these coordinates the background Minkowski metric takes the form 



/ 1 \ 



— sinh^y 

V — sinh^y sin^ ■»? / 



(27) 



and the velocity componets become = (1,0,0,0). Hence, the new coordinate frame is 
comoving. The metric (27) corresponds to the Milne cosmological model - a homogeneous, 
isotropic, expanding universe with the cosmological scale a — t and negative spatial curva- 
ture. 

The functional dependence of T on r follows from the energy- momentum conservation. 
For a perfect relativistic fluid the energy- momentum tensor is given by 

T/uu = (p + p)ui^u^ - pQu,^, (28) 

where p and p denote respectively the pressure and the energy density of the fluid. Prom 
the energy-momentum conservation 

T^""".^ = (29) 

applied to (28) we flnd 

^^V,m + (p + pK;m = 0, (30) 

where the subscript ; fj, denotes the covariant differentiation associated with the background 
metric. Since our fluid is dominated by massles pions at nonzero temperature, it is a rea- 
sonable approximation to assume the equation of state p = p/3 of an ideal gas of massless 
bosons. Then, Eq. (30) in comoving coordinates reads 

^ + ^ = (31) 

OT T 



with the solution 

P = Po(-) . (32) 



T . 

This expression combined with the density of the pion gas [31] 



P = (33) 

implies the temperature proflle 

T = To — . (34) 

T 

The constants Tq and Tq may be flxed from the phenomenology of high energy collisions. 
For example, if we choose Tq = IGeV, then a typical value of p = IGeV/fm^ at r ~ 5fm 
[32] is obtained with tq = 1.5 fm. In our case, with these values the interesting range of 
temperatures T between 100 and 200 MeV corresponds to r between 15 and 30 fm. In the 
following we work with Tq = IGeV and keep tq unspecifled so that physical quantities of 
dimension of time or length are expressed in units of tq. 



8 



3 Analog cosmology 



In this section we turn to study the analog metric and formation and properties of the 
apparent horizon in an expanding chiral fluid. To this end we outline the formalism in the 
first subsection and derive a condition for the apparent horizon for a general hyperbolic 
spacetime. In Sec. 3.2 we derive the analog metric for the expanding chiral fluid and study 
the properties of the analog apparent horizon. Then, in Sec. 3.3 we exploit the Kodama- 
Hayward definition of surface gravity to derive the Hawking temperature as a function of 
the parameters of the chiral fiuid, in particular, as a function of the local fiuid temperature. 



3.1 Radial null geodesies 

To study the apparent horizon in an expanding chiral fiuid we need to examine the behavior 
of radial null geodesies of the analog metric which we shall derive in Sec. 3.2. With hindsight, 
we first consider a space time of the form 

(is2 = p^rfdr'^ - a{Tf{dy^ + sinh^y dQ^), (35) 

where (3 and a are arbitrary functions of r. The metric tensor is 



— ct^ sinh^y 

V —a^ sinh^y sin^ i} J 



(36) 



This metric represents the class of hyperbolic (A; = — 1) FRW spacctimcs including the flat 
spacetime example (27). We denote by l'^ and the vectors tangent to outgoing and ingoing 
affinely parameterized radial null geodesies normal to a sphericall two-dimensional surface 
S. The tangent vectors are null with respect to the metric (36), i.e., 

G^^Al'^ = G^J'll'i = 0. (37) 

Using the geodesic equation 

l^V^l" = 0, (38) 

where the symbol V^^^ denotes a covariant derivative associated with the metric (36), one 
easily finds the tangent null vectors corresponding to four types of radial null geodesies, 

li^^q^a-^ (^p-\±a-\0,0) , (39) 

tangent to future directed and 

11 = q±a-^ ±«-\ 0, O) , (40) 

to past directed null geodesies , where g+, g_, q+, and g_ are arbitrary positive constants. 
The corresponding affine parameters A+ and A_ for the outgoing and ingoing null geodesies, 
respectively, are found to satisfy 

^ = -aP (41) 
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for future directed and 

dr q± 

for past directed null geodesies. For simplicity, from now on we set g+ = g_ = g+ = g_ = 1. 

The null vectors Z+ and l'^ point towards increasing and decreasing y, respectively. Hence, 
we adopt the usual convention and refer to l'^ (It) and the corresponding null geodesies as 
outgoing (ingoing) although increasing (decreasing) y does not necessarily imply increasing 
(decreasing) of the radial coordinate r. As we move along a geodesic the changes of the 
coordinates r and y are subject to the condition = of radial null geodesies, i.e., 

Oi 

dr = ±-dy (43) 

along the geodesic. Here the signs determine whether y is increasing or decreasing as we 
move along the geodesic. For example, for future directed null geodesies, it follows from (41) 
and (43) that an outgoing geodesic is directed along increasing y, i.e., y increases with A+, 
whereas an ingoing geodesic is directed along decreasing y, i.e., y decreases with A_. 

The key roles in the study of trapped surfaces are played by the expansion parameters 
£+ and £_ 

e± = V^l^ (44) 

of outgoing and ingoing null geodesies, respectively. Particularly important are the values 
of £+ and £_ and their Lie derivatives 

|;-fta,.- (45) 

in the neighborhood of a marginally trapped surface. As we shall shortly demonstrate, the 
relevant marginally trapped surface in the expanding chiral fluid is future inner marginally 
trapped. According to our convention described in Appendix A, a two-dimensional surface 
H is said to be future inner marginally trapped if the future directed null expansions on H 
satisfy the conditions: = 0, I'td^Sj^ln > and < 0. The future inner marginally 

trapped surface is the inner boundary of a future trapped region consisting of trapped 
surfaces with negative ingoing and outgoing null expansions. From now on we refer to this 
surface as the apparent horizon. 



From (39) and (40) we find 



^^--Al5±7-r7. (46) 




for future directed and 

for past directed radial null geodesies, where the overdot denotes a partial derivative with 
respect to r. The respective Lie derivatives are given by 




1 



tanh y J oP- 



2a 

e±: (48) 
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for future directed and 



a ajS \ tanh^i/ J o? 



for past directed radial null geodesies. 

For a spherically symmetric spacetime, the condition that one of the null expansions 
vanishes on the apparent horizon H is equivalent to the condition that the vector n^, normal 
to the surface of spherical symmetry, is null on H. In other words, the condition 

Vi,F\h^O, (50) 
where denotes either I'j^ or I'j^, is equivalent to the condition 

G^X^.k = 0. (51) 
This may be seen as follows. For the metric (36) the normal is given by 

= d^{asmhy). (52) 
The expansion £+ (or £_) defined in (44) may be written as 

V^P = ^i=a^(^=GF) = ^a^(V^P) + -^—Fn^ (53) 



V— G \J^^i ^ CKsinhy 

where h denotes the determinant of the metric 

V=(C (54) 

of the two-dimensional space normal to the surface of spherical symmetry. It may be shown 
that the first term on the right hand side of (53) vanishes identically by the geodesic equation. 
Hence, the expansion V^/'^ vanishes on H if and only if 

i%k = 0. (55) 

Suppose one of the expansions vanishes on if, i.e., Eq. (55) holds for either l\ or Since 
is null and both and are normal to B. and hence tangent to the two-dimensional space 
(r, y) with the metric (54), Eq. (55) implies hapn°'n^\H — 0. Hence, V fjl'^\H — implies 
G^'n^n^li? = 0. 

To prove the reverse it is sufficient to show that Z+ri/j 7^ and l^^n/i 7^ implies hai3n°'n^ 7^ 
0, which may be easily shown for a general two-dimensional metric in diagonal gauge. Then, 
the following statement holds: the vanishing of hapTi^iiP on H implies either /^n^|j^ = or 
l-np\H = 0. This together with (53), implies either £+|ff = or £_|ij = 0. 

Either from (50) or from (51) one finds the condition for the apparent horizon 

? ± = (56) 

(3 tanhy ^ ^ 
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3.2 Analog horizons 



Next we derive the analog metric and define the analog Hubble and the apparent horizons. 
Equation (12) may be written in the form [14, 33, 34] 



1 



_ CnduTT + -^Via, 7r)7r = 0, 

with the analog metric tensor, its inverse, and its determinant given by 



G 



HI/ 



-[Qh,/- (1 -4)m^m^], 



(57) 



(58) 



G 



C7I- 

a 



g^^ - (1 - -)u^vr 



(59) 



(60) 



Hence, the pion field propagates in a (3+1) -dimensional effective geometry described by the 

metric G^^. 

In the comoving coordinate frame defined by the coordinate transformation (26) the 
velocity is = (1, 0, 0, 0) and, as a consequence, the analog metric (58) is diagonal 



— r 



-T^ sinh^T/ 



V 



-T^ sinh^y sin^ 1? / 



(61) 



Here, the parameters a and c,^ are functions of the temperature T which in turn is a function 
of r. In the following we assume that these functions are positive. The metric (61) is precisely 
of the form (36) with 



'ac-n 




(62) 



The physical range of r is fixed by Eq. (34) since the available temperature ranges between 
T = and T = T^. Hence, the proper time range is Tc < r < 00 where the critical value Tc 
is related to the critical temperature as Tc/tq = Tq/Tc. The metric is singular at r = Tc. 

In Fig. 3 we plot the expansions and S- of outgoing and ingoing radial null geodesies, 
respectively, as functions of r for an arbitrarily chosen fixed time t — 6to and, similarly, 
as functions of y for a fixed r = 5.77tq. In the lower two panels we plot the derivative of 
the outgoing null expansion e+ along the ingoing null geodesic. The outgoing null expansion 
decreases with increasing r from positive to negative values and vanishes at the point r = th, 
whereas the ingoing null expansion remains negative. At this point the derivative of with 
respect to A_ is positive. According to the standard convention described in Appendix A 
the region {r > rn, t — 6to} is future trapped and the location rn marks its inner boundary. 
Thus, the sphere at Th is future inner marginally trapped. 

Spacetimc diagrams corresponding to the metric (61) are presented in Fig. 4 showing 
future directed radial null geodesies. The origin in the plots in both panels corresponds to 
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Figure 3: Null expansions e+ and S- as functions of r (top left panel) and y (top right panel) 
for fixed t = 6ro and fixed r = 5.77ro, respectively. Similarly, the bottom panels depict the 
derivative of as functions of r and y for fixed t and r, respectively. 



the critical value Tc at which c-^ vanishes. Numerically, with the chosen Tq = 1 GeV we have 
Tc/tq = 5.47. The geodesic lines are constructed using 

y = ± dT'cT,{T')/T' + const (63) 



which follows from (43) with (62). As mentioned in Sec. 3.1, increasing (decreasing) y does 
not necessarily imply increasing (decreasing) of the radial coordinate r. With the help of 
the coordinate transformation (26) the shift dy along a geodesic may be expressed in terms 
of dr 

dy = 7 — 7— dr. (64) 

where we have used (25) and (62). We note that if f < c^r, increasing(decreasing) y corre- 
spond to increasing (decreasing) r for both signs, whereas if v > c-„ increasing y corresponds 
to increasing r for an outgoing and decreasing r for an ingoing geodesic. 
Using (63) we introduce null coordinates 

^ = l{y + £ dr'P{T')/a{r')^ ; « = ^ + [ dr' P{T')/a{r')^ , (65) 
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5.5 6.0 6.5 

Figure 4: Spacetime diagram of outgoing (full line) and ingoing (dashed line) radial null 
geodesies in (r, y) coordinates. The shaded area represents the evolution of the trapped 
region. The trapping horizon is represented by the full bold line with the endpoint at 
T = Tmax = 6.0182ro. The dashed and dash-dotted bold lines represent the evolution of the 
analog and naive Hubble horizons, respectively. 



ranging in the intervals [0, +oo) and (— oo, +oo), respectively, with a condition < w ± u < 
oo. In these coordinates the metric (35) becomes 

ds^ = (Adudw - sinh^ - u) dQ"^) . (66) 

The singularity at r = Tc is mapped onto the entire u + w = line. Next, we compactify the 
spacetime using the coordinate transformation 

W = ta.nh.w] f/ = tanhM. (67) 

The coordinates W and U range in the intervals [0,1] and [—1,1], respectively, with a 
condition < W ± U < 2. Furthermore, the rotation 

T = ^{W + Uy, R = 1=(W-U), (68) 



brings the metric to a conformally flat form 

ds' 



;i-f/2)(i-ty2) 



dT^ - dR^ - R^ dVL' 



(69) 
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Figure 5: Conformal diagram for the analog spacetime described by the metric (61). The 
Unes parallel to the W and U axes, are future directed outgoing and ingoing null geodesies, 
respectively. The singularity at r = Tc is represented by the wavy line and the apparent 
horizon by the full line. Te shaded area in between represents the evolution of the trapped 
region. The spacelike and future timelike infinities are denoted by io and i+, respectively. 
The future null infinity X+ is represented by the line T + R = 



where both coordinates R and T range in the interval [0, v^] with a condition R + T < 
a/2. The conformal diagram representing our analog spacetime is depicted in Fig. 5. The 
singularity at r = Tc is mapped onto the segment [0, a/2] on the horizontal axis. 
The coordinate transformation 

t' = ( l3dT (70) 



brings the metric (61) to the standard form of an open k = —1 FRW spacetime metric with 
the cosmological time t'. The time coordinate t' is related to the original time t via r and the 
transformation (26). The analog cosmological scale is a{t') = a(r(t'))/ro, where the costant 
ro is related to the spatial Gaussian curvature K = — I/tq. The proper distance is dp = ay 
and the analog Hubble constant 

ap 

where the overdot denotes a partial derivative with respect to r. Then, we define the analog 
Hubble horizon as a two-dimensional spherical surface at which the magnitude of the analog 
recession velocity 

Wrec = "Hrfp = (72) 



15 



equals the velocity of light. Hence, the condition 



defines the location of the analog Hubble horizon. Note that the radial fluid velocity v in (25) 
and the analog recession velocity (72) are quite distinct quantities - in an expanding fluid v 
is allways positive and less than the velocity of light c = 1 whereas free may be positive or 
negative depending on the sign of a and its magnitude may be arbitrary large. 

A two-dimensional spherical surface on which the radial velocity v equals the velocity of 
pions deflnes another horizon, which we refer to as the naive Hubble horizon. This horizon 
is obviously distinct from the analog Hubble horizon deflned above. The evolution of the 
naive and the analog Hubble horizons with r arc depicted in Fig. 4. 

Next we introduce the concept of analog marginally trapped surface or analog apparent 
horizon following closely the general considerations of Sec. 3.1 and Appendix A. Formation 
of an analog apparent horizon in an expanding hadronic fluid is similar to the formation 
of a black hole in a gravitational collapse although the role of an outer trapped surface is 
exchanged with that of an inner trapped surface. Unlike a black hole in general relativity, 
the formation of which is indicated by the existence of an outer marginally trapped surface, 
the formation of an analog black (or white) hole in an expanding fluid is indicated by the 
existence of a future or past inner marginally trapped surface. 

Equation (56) with (62) defines a hypersurface which we refer to as the analog trapping 
horizon. Any solution to Eq. (56) , e.g., in terms of r for fixed t, gives the location of the 
analog apparent horizon rn- For example, the radius rn = 1.53ro computed using (56) for 
fixed t — 6to is the point of vanishing outgoing null expansion which marks the location of 
the apparent horizon in the top left panel of Fig. 3. From (46) it follows that the region of 
spacetime for which 



is trapped. It is future trapped if d < and past trapped if a > 0. The condition (74) 
can be met only if \/3/a\ < 1 which holds for r between Tc and Tmax- At r = Tmax we have 
1/3/ a\ = 1 so the endpoint of the trapping horizon in Fig. 4 is at r = Tmax, tanhy = 1. 

We find that a is negative for r in the entire range Tc < r < Tmax and, according to 
(71), the analog Hubble constant is always negative. Hence, our analog cosmological model 
is a contracting FRW spacetime with a negative spatial curvature. The shaded area left 
of the bold line in Fig. 4 represents the time evolution of the future trapped region. Note 
that the analog Hubble horizon is always behind the apparent horizon whereas the naive 
Hubble horizon may be located ahead of or behind the apparent horizon depending on 
the magnitude of a. The naive Hubble and apparent horizons coincide if a and are r 
independent constants. 

The apparent horizon is generally not a Killing horizon and normally does not coincide 
with the event horizon (one exception is de Sitter spacetime). Moreover, the apparent 
horizon exists in all FRW universes [35] , whereas the event horizon does not exist in eternally 
expanding FRW universes with the equation of state w > —1/3 (see, e.g., [36]). For the 
metric (61), the event horizon is defined by 



tanhy > 1(3/ a 



(74) 




(75) 
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In our chiral fluid model the integral on the righthand side diverges at the upper limit because 
— > 1 as T ^ oo and hence, the analog event horizon does not exist. In contrast, as we 
have demonstrated, the analog apparent horizon does exist. 



3.3 Analog Hawking effect 

One immediate efi^ect related to the apparent horizon is the Hawking radiation. Unfortu- 
nately, in a non-stationary spacetime, the surface gravity associated to the apparent horizon 
is not uniquely defined [37]. Several ideas have been put forward how to generalize the defi- 
nition of surface gravity for the case when the apparent horizon does not coincide with the 
event horizon [38, 39, 40, 41]. In this paper we use the prescription of [41] which we have 
adapted to analog gravity in our previous paper [16]. This prescription involves the so called 
Kodama vector [42] which generalizes the concept of the time translation Killing vector 
to non-stationary spacetimes. The analog surface gravity n is defined by 

K^l^=daiV^h'^^kn0), (76) 



where the quantities on the righthand side should be evaluated on the trapping horizon. The 
metric ha/3 of the two-dimensional space normal to the surface of spherical symmetry and 
the vector normal to that surface are given by (54) and (52), respectively. 

The definition (76) differs from the original expression for the dynamical surface gravity 
[41, 43] by a normalization factor k which we have introduced in order to meet the require- 
ment that should coincide with the time translation Killing vector for a stationary 
geometry. For the metric (61) with (62) we have found [16] 

k — j3 ^cosh^ y — sinh^ y — ^ . (77) 



Then, the definition (76) yields 



K 



^ a{a^ + aa- ^^) 2^^ (rd -a)v + {ad? - Ito? /3'^ra)v'^ (78) 



evaluated on the trapping horizon. The above expression may be somewhat simplified by 
making use of the horizon condition (56). We find 

^ 1 + 2c^-t;(l -v)-{2 + c^)v^ a v 

~ 2r -fv{l + c^vf 2/3 7(1 + c^v)^ ^ ' 

evaluated on the trapping horizon. 

It is worthwhile analysing the limitting case of (78) when the quantities a, and are 
constants. Then d = a/r, a = {) and the apparent horizon is fixed by the condition v = Ct^. 
At any chosen time t = t(1 — c^)~^/^ the horizon is located at rn = c^^t and the expression 
for K, reduces to 

«=1 ^ (80) 
2t 2t ^ ^ 
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Hence, the analog surface gravity is finite for any physical value of and is maximal when 
= 0. However, with = the horizon degenerates to a point located at the origin r = 0. 
In the left panel of Fig. 6 we plot k as a function of r as given by (79). The corresponding 
temperature defined as 

r„ = - (81) 

represents the analog Hawking temperature of thermal pions emitted at the apparent horizon 
as measured by an observer near infinity. Since the background geometry is flat, this tem- 
perature equals the locally measured Hawking temperature at the horizon. Thus, equation 
(81) with (78) corresponds to the fiat spacetime Unruh effect. 

As we move along the trapping horizon the radius of the apparent horizon increases and 
the Hawking temperature decreases rapidly with r. Hence, there is a correlation between 
Th and the local fluid temperature T which is related to r by (34). In the right panel of 
Fig. 6 we show the Hawking temperature Th as a function of the fluid temperature T at the 
apparent horizon. 

In our previous paper [16] we have shown that the surface gravity diverges as 

n = {r] + l/2){r-T,)~' (82) 

at the singular point, where 77 is a constant related to the scaling of the quantity ^J'a/c^ 

oc (Te-T)-" (83) 

in the neighborhood of the critical point. The constant rj may be roughly estimated as 
follows. The estimate of the function S(g) defined in (19) in the neighborhood of = for 
small a yields [24] 

S(0,q^)~-g2; S(go,0)~^go' (§4) 
a a 
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By comparing this with (18) we deduce the behavior of the quantities a and b for small a 




a/c^ diverges as a/cjr oc cr From Eq. (8) we find a oc (Tc — T)^'^ near the critical point 
which yields r] = 3/8. 



Numerically, by fitting ya/cT^ in the close neighborhood of Tc to the function (83) with 
the critical temperature Tc — 182.822 MeV obtained numerically from (8), we find rj — 0.253. 
A more refined analysis based on scaling and universality arguments of Son and Stephanov 
[23] yields rj = 0.1975 [16]. 



We have demonstrated that, owing to the analog gravity effects in high energy collisions, a 
close analogy may be drawn between the evolution of a hadronic fiuid and the spacetime 
expansion. Using the formalism of relativistic acoustic geometry we have analyzed the ex- 
panding chiral fluid in the regime of broken chiral symmetry. The expansion which takes 
place after the collision is modelled by spherically symmetric Bjorken type expansion. The 
propagation of massless pions in the chiral fluid provides a geometric analog of expanding 
spacetime equivalent to an open (k — —1) FRW cosmology. The geometry depends on the 
parameters a and b of the effective Lagrangian defined in Sec. 2. The elements of the analog 
metric tensor are functions of the spacetime coordinates via the temperature dependence of 
a and the pion velocity c^. The pions propagate slower than light with a velocity close to 
zero in the neighborhood of the critical point of the chiral phase transition. 

A trapped region forms for radial velocities of the fiuid beyond the value defined by 
Eq. (56). This value defines a hypersurface shown in Fig. 4 which we refer to as the analog 
trapping horizon, at which the outgoing radial null expansion vanishes. Our trapping horizon 
is foliated by future inner marginally trapped surfaces and is equivalent to the trapping 
horizon in contracting FRW spacetime, i.e., in dynamical spacetime with a negative Hubble 
constant. The shaded area in Fig. 4 represents the time evolution of the future trapped 
region, with the future inner marginally trapped surface (or the future apparent horizons) 
as its inner boundary. This marginally trapped surface may be regarded as an "outer" white 
hole: the ingoing pions (future directed ingoing null geodesies) freely cross the apparent 
horizon whereas the outgoing cannot penetrate the apparent horizon. This is opposite to an 
expanding FRW universe where the inner marginally trapped surface acts as a black hole: 
the future directed ingoing null geodesies cannot escape the apparent horizon whereas the 
outgoing null geodesies freely cross the apparent horizon. 

We have studied the Hawking effect associated with the analog apparent horizon using 
the Kodama- Hay ward definition of surface gravity adapted to the analog gravity geometry. 
The Hawking temperature correlates with the local temperature of the fiuid at the appar- 
ent horizon and diverges at the critical point. In contrast to the usual general relativistic 
Hawking effect, where the Hawking temperature is tiny compared with the temperature of 




4 Summary and discussion 
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the background, the analog horizon temperature is of the order or even larger than the local 
temperature of the fluid. 

The analog Hawking radiation of pions should not be confused with the Hawking-Unruh 
radiation of hadrons of Castorina ct al. [44]. The latter is a usual Unruh effect due to the 
acceleration of quark-antiquark pairs produced in particle collisions, whereas the former is 
an analog thermal radiation due to effective geometry of the chiral fluid. 

A spherically symmetric Bjorken expansion model considered here may be phcnomcno- 
logically viable as a model of hadron production in e+e^ but it is certainly not the best 
model for description of high energy heavy ion collisions. It would be desirable to apply 
our formalism to a more realistic hydrodynamic model that involves a transverse expansion 
superimposed on a longitudinal boost invariant expansion. In this case the calculations be- 
come rather involved as the formalism for general nonspherical spacetimes is not yet fully 
diveloped. This work is in progress. 

In conclusion, we believe that the study of analog gravity in high energy collision may in 
general improve our understanding of both particle physics phenomenology and dynamical 
general relativistic systems. 
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A Trapped surfaces in general relativity 

Following Hay ward [40] we summarize here the relevant definitions related to trapped sur- 
faces. 

• Trapped surface. Let E denote a spacelike hypersurface, e.g., a hypersurface defined 
by t — const. A two-dimensional surface 5" C E with spherical topology is called a 

trapped surface on E if the families of ingoing and outgoing null geodesies normal to 
the surface are both converging or both diverging. More precisely, the null expansions 
£+ = I'^-n and £_ = on a trapped surface S should satisfy > 0. One 

distinguishes between a past trapped surface for which both £+ and £_ are positive, 
and a future trapped surface for which both and £_ are negative. 

• Trapped region. A set of future trapped surfaces (or closed trapped surfaces [35, 45]) 
on E is referred to as a future trapped region. Similarly, a set of past trapped surfaces 
on E is called a past trapped region. 

• Mcirginally trapped surface. A two-dimensional surface H is said to be marginally 
trapped if one of the null expansions vanishes on H, i.e., if either s+\h = or e^ln = 0. 
A marginally trapped surface is also referred to as an apparent horizon although, strictly 
speaking, the original definition of Ellis and Hawking [45] involves in addition the 
assumption of asymptotic flatness. 
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• Outer mcirginally trapped surface. A surface H is said to be future (past) outer 
marginally trapped if on H the future (past) directed outgoing null expansion vanishes, 
its derivative along the ingoing null geodesic is negative, and the ingoing null expansion 
is negative, i.e., if e^ln = 0, < 0, and S-\h < 0. Equivalently, the future 
(past) outer marginally trapped surface may be defined as the surface on which the 
past (future) directed ingoing null expansion vanishes, its derivative along the outgoing 
null geodesic is negative, and the outgoing null expansion is positive, i.e., if £_|if = 0, 
l+dfj,e^\H < 0, and £+|// > 0. 

• Inner mcirginally trapped surface. A surface H is said to be future (past) inner 
marginally trapped if on H the future (past) directed outgoing null expansion vanishes, 
its derivative along the ingoing null geodesic is positive, and the ingoing null expansion 
is negative, i.e., if = 0, l'td^e+\H > 0, and E-Ih < 0. Equivalently, the future 
(past) inner marginally trapped surface may be defined as the surface on which the 
past (future) directed ingoing null expansion vanishes, its derivative along the outgoing 
null geodesic is positive, and the outgoing null expansion is positive, i.e., {s-\h — 0, 

> 0, and e+Ih > 0). 

• Trapping horizon. A hypersurface foliated by inner or outer marginally trapped 
surfaces is referred to as an inner or outer trapping horizon, respectively. 

According to this classification we distinguish four physically relevant classes: 

1. Future outer marginally trapped surface is the outer boundary of a future trapped 
region and is typical of a black hole. 

2. Past outer marginally trapped surface is the outer boundary of a past trapped region 
and is typical of a white hole. 

3. Future inner marginally trapped surface is the inner boundary of a future trapped 
region and represents an "outer" white hole. This situation is physically relevant in 
the cosmological context for a shrinking universe, i.e., for an FRW spacetime with 
d < 0. 

4. Past inner marginally trapped surface is the inner boundary of a past trapped region 
and represents an "outer" black hole.This situation is physically relevant in the context 
of an expanding FRW universe [35, 46]. 
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